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Based on a numerical ab initio study, we discuss a structure model for a broad boron sheet, 
which is the analog of a single graphite sheet, and the precursor of boron nanotubes. The sheet has 
linear chains of sp hybridized a bonds lying only along its armchair direction, a high stiffness, and 
anisotropic bonds properties. The puckering of the sheet is explained as a mechanism to stabilize the 
sp a bonds. The anisotropic bond properties of the boron sheet lead to a two-dimensional reference 
lattice structure, which is rectangular rather than triangular. As a consequence the chiral angles of 
related boron nanotubes range from 0° to 90°. Given the electronic properties of the boron sheets, 
we demonstrate that all of the related boron nanotubes are metallic, irrespective of their radius 
and chiral angle, and we also postulate the existence of helical currents in ideal chiral nanotubes. 
Furthermore, we show that the strain energy of boron nanotubes will depend on their radii, as well 
as on their chiral angles. This is a rather unique property among nanotubular systems, and it could 
be the basis of a different type of structure control within nanotechnology. 



PACS numbers: 81.07.De, 73.63.Fg, 61.46,+w 

I. INTRODUCTION 

Boron is an electron deficient element which has a 
rather fascinating chemistry. Pure boron compounds nei- 
ther have a purely covalent nor a purely metallic charac- 
ter. This results in a chemical versatility, which is unique 
among the elements of the periodic table. 

The classical bulk modifications of boron are based on 
B12 icosahedra. The simplest boron phase is rhombohe- 
dral a-boron, 2 where boron icosahedra are centered on 
the edges of a rhombohedral unit cell. A different pic- 
ture arises for boron clusters, where quasiplanar isomers 
turn out to be more stable than their icosahedral coun- 
terparts. This is the consequence of an Aufbau princi- 
ple for elemental boron clusters postulated by Boustani. 3 
This Aufbau principle generally states that stable boron 
clusters can be constructed from two basic units only: a 
pentagonal pyramidal B6 unit and a hexagonal pyramidal 
B7 unit, 3 and it implies quasiplanar, 4 tubular, 5,6 convex, 
and sphericafi boron clusters. The existence of quasi- 
planar clusters or "sheets" was recently confirmed by 
experiment, 8 in perfect agreement with earlier theoret- 
ical predictions. 9 Furthermore, the existence of quasipla- 
nar boron clusters implies the formation of boron nano- 
tubes and/or boron fullerenes, because during synthesis, 
a growing (quasi-) planar boron cluster tends to remove 
dangling bonds by forming closed tubular or polyhedral 
modifications. And indeed, recent experimental studies 
demonstrate the existence of boron nanotubes. 10,11 

Carbon nanotubesi^ on the other hand are a structural 
paradigm for all nanotubular materials and they can be 
seen as cylindrical modifications of graphite, which may 
geometrically be constructed by cutting a rectangular 
piece out of a single graphene sheet and rolling it up to 



form a tube. Almost all properties of carbon nanotubes 
can be derived from the properties of a single graphene 
sheet, which means that a profound understanding of 
graphite is the key to understand the basic properties 
of carbon nanotubes. The same relation holds for boron 
sheets (BSs) and boron nanotubes (BNTs): understand- 
ing the structure and the properties of BSs will be crucial 
for our understanding of the basic properties of BNTs. 

This paper builds on previous workiiii to establish 
such a basic connection between BSs and BNTs, but it 
should be pointed out that our previous reasoning was 
mainly based on the individual structures of finite sized 
quasiplanar boron clusters Ai 5 ^^ Using ab initio struc- 
tural optimization methods for solid systems we could 
finally discriminate among different structure models for 
layered boron compounds and establish a simple model 
for a broad and stable BS. 

After a detailed description of this search process, we 
will analyze the properties of the most stable structure 
model. Then we will show how these results may be used 
to explain the structure, the stability, the electronic and 
the mechanical properties of BNTs. In particular the 
somewhat surprising constriction of zigzag BNTs, which 
has been reported in a recent publication^ may now be 
clearly understood on the basis of the elastic properties 
of BSs. 

It must be pointed out that up to now, a broad BS, 
which would be the analog of a single graphite sheet, 
could not be found experimentally. But when writing up 
this paper we became aware of an interesting work by 
Evans et al^ who consider three BS models and five 
BNTs of small tube radii, and the work of Cabria et alm^ 
who study two BS models and three BNTs. Although our 
results are certainly based on a more extensive search for 
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stable BS and BNTs, our findings for the stable BS are, 
from a structural and energetic point of view, in excellent 
agreement with these authors. Thus the present struc- 
ture model could independently be confirmed by three 
different groups. However, there is still some disagree- 
ment about the ground state structures of BNTs. Lau 
et a/., 20 for example, have recently reported about struc- 
tures for BS and BNTs, which are very different from the 
structure models of Evans et al and Cabria et a/., but 
the present study is in clear favor of the latter. 

II. METHODS 

As pointed out by Pauling 1 elemental boron has a 
complicated and rather versatile chemistry. Therefore 
the only reliable theoretical tools, which may allow for a 
proper description of boron chemistry, are first principles 
calculations. 3 

In order to carry out structural optimizations of 
BSs and BNTs we used the VASP package, version 
4. 4. 6. 21-22 The latter is a density functional theory 23 
based ah initio code using plane wave basis sets and 
a supercell approach to model solid materials, sur- 
faces, or clusters™ During all of our simulations, the 
electronic correlations were treated within the local- 
density approximation (LDA) using the Perdew-Zunger- 
Ceperley- Alder exchange-correlation functional and 
the ionic cores of the system were represented by ul- 
trasoft pseudopotentials 27 as supplied by Kresse and 
Hafner™ The /c-space integrations were carried out us- 
ing the method of Methfessel and Paxton 2 ^ in first order, 
where we employed a smearing width of 0.3 eV. 

With the help of the VASP program, one can deter- 
mine interatomic forces, which may be used to relax the 
different degrees of freedom for a given decorated unit 
cell. Eventually one will detect some atomic configura- 
tions, which correspond to (local) minima on the total 
energy landscape. In order to carry out those extensive 
structure optimizations in a more effective way, we em- 
ployed a conjugate gradient algorithm, 24 and we allowed 
all of the atomic coordinates to relax, as well as all but 
one lattice parameter. This rigid lattice parameter would 
fix the inter layer separation for BS and the intertubular 
distance for BNTs at 6.4 A, which effectively makes them 
stand-alone objects. The sizes of the fc-point meshes for 
different systems with different unit cells were individ- 
ually converged, such that changes in the total energy 
were reduced to less than 3 meV/atom. In the course 
of a structural optimization run, all interatomic forces 
were finally reduced to less than 0.04 eV/A. The cutoff 
energy for the expansion of the electronic wave functions 
in terms of plane waves was 257.1 eV for the relaxation 
runs, and 321.4 eV for a final static calculation of the 
total energy. 

The cohesive energies given in Tables [I] and [0] were 
calculated from 

i^coh = Ebi n( ±/n. (1) 




FIG. 1: (Color online) Top view of a quasiplanar boron sheet. 
In a planar projection the atoms form an almost perfect trian- 
gular lattice. The basic structural unit is a hexagonal pyrami- 
dal B7 cluster, as suggested by the Aufbau principle (Ref. 0) 
(see text). 



^bind is the the atomic binding energy per unit cell and 
n is the number of atoms per unit cell. Therefore in our 
definition £? co h will be a positive number. 

For band structures and the analysis of Fermi sur- 
faces in Sec. 1111 B 21 and Appendix [B] we used the 
Stuttgart TB-LMTO-ASA package, which is a density 
functional theor^ 3 - based code using short range 30 lin- 
earized muffin-tin orbitals 3 ^ within the atomic sphere ap- 
proximation (ASA). It allows static calculations of the 
electronic properties for periodic systems. We used the 
non-spin polarized LDA exchange-correlation functional 
of Barth and Hedin 32 and a /c-mesh of 30 x 30 x 3. 



III. BORON SHEETS 

A. Finding a structure model 

Following the Aufbau principle 3 - a BS is basically 
a quasiplanar arrangement of hexagonal pyramidal B7 
units. A planar projection of such a system will always 
form some kind of triangular lattice (see Fig. QJ. How- 
ever, the out of plane modulation (i.e., the puckering) 
remains unspecified by the Aufbau principle. The latter 
has to be determined using ah initio structural optimiza- 
tions, after setting up a suitable supercell that will allow 
for a systematic generation of various periodic puckering 
schemes. 

The versatile chemistry of boron is reflected in a com- 
plicated energy landscape, which is full of local min- 
ima. Therefore the standard optimization techniques like 
the conjugate gradients method used in this study are 
most likely to find local minima, rather than global min- 
ima. Therefore we examined the energy landscape quite 
carefully by performing many optimization runs, which 
started from quite diverse initial configurations. 

The basic puckering schemes were taken from the 
structures of B22, B32, and B46 clusters, which are de- 
scribed in Ref. 4. We repeated the puckering periodi- 
cally in a triangular supercell containing 16 atoms (see 
Figs. Efb)-|2Id)), and optimized the resulting structures. 
For the sake of comparison we also examined a fiat BS 



3 





FIG. 2: (Color online) Different structure models for broad 
boron sheets. Each supercell (thin lines) contains 16 atoms, 
(a) A simple flat sheet is metastable. (b) A simple up and 
down puckering seems to be the most stable modulation. 
Structures (c) and (d) are unstable. Models (b), (c), and 
(d) are periodic repetitions of structural motives taken from 
B22, B46, and B32 clusters, described in Ref. 0. 



(see Fig. |2Ja)). 44 The flat boron sheet (a) occupies a lo- 
cal minimum on the energy landscape with a cohesive 
energy of 6.76 eV/atom, but small out-of-plane elonga- 
tions of individual atoms immediately cause a puckering 
of the BS. This was confirmed by shifting one atom 0.1, 
0.2, and 0.4 A out of plane and reoptimizing the resulting 
structures. Thus model (a) turns out to be metastable 
(as also pointed out in Refs. and EH) : any thermal vi- 
bration would lead to a permanent deformation of a flat 
boron sheet. Models (c) and (d) are completely unsta- 
ble, and they both relax to structure (b). In order to 
scan the energy landscape for other candidate structures 
we took sheet (a) and shifted each of the 16 atoms out of 
the plane, employing a random elongation Az between 
+0.4 and —0.4 A. Those structures were reoptimized as 
before. It turns out that 8 out of 11 optimizations led 
to model (b), while the remaining three runs resulted in 
a metastable kinked structure with a cohesive energy of 
6.86 eV/atom (see Appendix [HJ). 

The fact that models (c) and (d) as well as 8 out of 
11 randomly puckered sheets would relax to model (b) 
means that structure (b) defines a rather pronounced 
minimum on the energy landscape. The high structural 
stability of model (b) is confirmed by its high cohesive 
energy of 6.94 eV/atom, which is the highest cohesive en- 
ergy of all BSs that we found. We thus conclude that the 
most suitable structure model for a broad BS will be (b), 
being 0.18 eV/atom more stable (0.21 and 0.26 eV/atom 
in Refs. andUM respectively) than an unrealistic flat 
BS. The puckering itself seems to be an important mech- 
anism to stabilize the BS, 19 which will be examined in 
more detail in Sec. IIII Bl 



FIG. 3: (Color online) The orthorhombic unit cell of model 
(b) with two basis atoms (see Table UJ. In a ^-projection 
atom 1 is located at the corners of a rectangular unit cell, 
while atom 2 is located at the center of the unit cell. Along 
the z direction the boron atoms will generate a simple up and 
down puckering, with puckering heights around Az = 0.82 A. 



In order determine the lattice structures of (a) and (b) 
we performed LDA calculations, where we would fix the 
unit cell of each system for a series of Cartesian lattice 
constants A or B, whereas all of the internal (atomic) 
degrees of freedom were allowed to relax. The resulting 
total energies for a given set of lattice constants were fit- 
ted to polynomial curves E(A) and E(B), from which 
we determined the equilibrium properties of the systems. 
The results are summarized in Table [IJ The diagonal el- 



TABLE I: Detailed LDA description of the optimized lattice 
structures of the flat (a) and puckered (b) boron sheets (see 
Figs. |21 and EJ, their bond lengths, cohesive energies E co h 
(Eq. Q), and their elastic moduli C x = Cu and C y = C22 
obtained after stretching a sheet along the Cartesian x or y 
direction (Eqs. J2J and @). 



Sheet 


(a) 


Flat 


(b) 


D uckered 


Lattice type 


Triangular (2D) 


Orthorhombic (3D) 


Lattice param. (A) 


A = 


= 1.69 


A = 


2.82 








B 


1.60 








C = 


arbitrary 


Primitive vectors 


ai 


ai = 


= A(1,0,0) 




a 2 




a 2 = 


= 5(0,1,0) 








a 3 = 


= C(0,0,1) 


Atoms/unit cell 


1 




2 




Atomic pos. (A) 


#1 


= (0,0) 


#1 = 


= (0,0,0) 








R2 


= (§^,§£,0.82) 


Bond lengths (A) 


<zb- 


-b = 1.69 


aS_ 


b = 1-60 






diat 

a B J 


g,„al = 


E coh (eV) 


6.76 


6.94 




Elastic modulus 


C x 


= Cy = 0.75 


C x z 


= 0.42 


(TPa) 






Cy z 


= 0.87 
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ements of the elastic tensor G X = G\\ and C y = C22 may 
be interpreted as a first approximation to a macroscopic 
Young's modulus. They were calculated as follows: 



A | 


fd 2 E{A) 


Bh 1 


V dA 2 




( d 2 E(B) 


— 

Ah 


\ dB 2 



h is the height of the BS, and it was defined as h = 
Az + 2i? v dw; is the puckering height of the sheet and 
^vdw is the van der Waals radius 4£ Aq and Bq are the 
equilibrium lattice constants. 

The optimized planar model (a) seems to form a tri- 
angular lattice with one atom per unit cell and a sin- 
gle lattice constant A, which is in the range of a typi- 
cal boron-boron bond length A = <2b-b = 1.69 A. But 
within the accuracy of the given methods, we cannot re- 
ally decide whether the lattice structure is perfectly tri- 
angular or slightly less symmetric. Assuming perfect tri- 
angular symmetry the two elastic moduli C x and C y are 
equal, and they are surprisingly big: C y = C y w 750 
GPa. Which means that even if the flat BS is metastable 
compared to other model boron sheets, it nevertheless 
has an extraordinary high stiffness. In Appendix [B] we 
will analyze model (a) in more detail. 

In Fig. |3| we depicted the unit cell of model (b). It 
consists of two basis atoms, and its planar projection 
is almost triangular, but not quite so. It is common 
to describe such a system with a face centered rectan- 
gular unit cell with lattice constants A and B. For 
A/B = y/3 = 1.732 a planar projection of the system 
would be equivalent to a triangular system. In our case 
A/B = 1.76, which is a small, but noticeable departure 
from triangular symmetry. Due to a puckering height 
of Az = 0.82 A, such a system might best be described 
using a three-dimensional orthorhombic unit cell. The 
corresponding lattice parameters and bond lengths can 
be found in Table [fl 46 



B. Properties of the model boron sheet 

In this section we will analyze the properties of model 
(b), which turns out to be the most stable structure for 
broad BS. Therefore, whenever we write 'boron sheet' 
(BS) in the following, we will only refer to model (b). 

In order to compare the BS with a known boron struc- 
ture we also calculated the cohesive energy of the a 
boron, which turns out to be 7.51 eV/atom. This corre- 
sponds to an energetic difference of 0.57 eV/atom (0.58 
and 0.57 eV/atom in Refs.H^and U^. respectively), which 
is huge, but one has to take into account that we are com- 
paring a single boron sheet with a bulk reference struc- 
ture. 



zigzag 




o-bonds 



FIG. 4: (Color online) Orange (gray): charge density con- 
tours of the boron sheet (model (b)) at 0.9 e/A 3 . One observes 
parallel linear chains of sp hybridized a bonds lying along the 
armchair direction. 



1. Mechanical properties 

The elastic modulus of model (b) strongly depends on 
the stretching directions. In Table [I] we roughly find that 
C y « 2C X . How can one explain those rather obvious 
anisotropics? 

To this end, let us have a look at the charge density 
of the BS (see Fig. UJ). We clearly observe some par- 
allel linear chains of a bonds lying along the armchair 
direction. Their bond length is ag_ B = 1.60 A. At lower 
densities (p < 0.7 e/A 3 , not displayed) a largely homoge- 
neous distribution with a rather complex shape appears, 
which may be assigned to multicenter bonding typical for 
boron materials. An analysis of the electron localization 
function 33 (ELF) leads to similar results, such that we 
obtain the following preliminary picture of the bonding: 
on a first level the sheet is held together by homogeneous 
multicenter bonds, but on a second level there are strong 
(7 bonds lying only along the armchair direction. 

Due to the strong a bonds, any stretching of the BS 
along the armchair (= y) direction will be much harder 
than a similar stretching along its zigzag (= x) direc- 
tion, where only the slightly weaker multicenter bonds 
are involved. These results are quite different from the 
results of Evans et al., who conjecture that the a bonds 
are strong but soft. 18 But here we clearly observe that 
the a bonds are strong and stiff. However other basic 
findings of Evans et al. are in good agreement with our 
results for flat and puckered BSs. 

In general the elastic moduli involved are quite high; 
the stiffness of the a bonds along the armchair direction 
is comparable to the stiffness of a graphene sheet. Fur- 
thermore the broken triangular symmetry of the BS's 2D 
lattice structure, as mentioned in Sec. IIII Al is another 
immediate consequence of the anisotropic bond proper- 
ties. 

Evans et al. also found that BNTs of different chirali- 
ties have different stiffnesses^ This can be confirmed by 
our bonding picture, although our results suggest that 
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FIG. 5: The band structure of the model BS. The fatness of 
the bands indicates their sp character, and it shows that the 
a bonds in Fig. 2|must be of sp type. The Fermi energy Ep 
lies at E — 0, G is the Y point. 



zigzag BNTs should be somewhat stiffer than armchair 
BNTs, while Evans et al. noted the opposite (the arm- 
chair and zigzag direction are swapped in their and our 
treatment, see Appendix IA lj) . We thus conclude that 
the relation between the microscopic elastic modulus and 
the macroscopic Young's modulus must be rather com- 
plicated in the case of BS and BNTs. 




FIG. 6: (Color online) The two-dimensional Fermi surface of 
the boron sheet. It consists of two contours in red (black) and 
yellow (gray) , which correspond to the two bands crossing the 
Fermi energy in Fig.E2 



2. Electronic properties 

The two-dimensional band structure of the BS 
E BS (k Xj k y ) is plotted in Fig. [5] along lines of high sym- 
metry. The BS is metallic, as there are two bands cross- 
ing the Fermi energy, which is in perfect agreement with 
earlier studies of BSs.Aiii 

In order to find out about the hybridization of the a 
bonds, we plotted the corresponding amount of s and p y 
character indicated by the fatness of the bands £ We do 
not find individual dispersions of s or p bands, and the 
lowest lying bands show dispersions which share s and 
p y character. That means they are bands consisting of 
sp hybridized orbit als: 

\*Pa) = ^=(|s) + |ft,» 

\m) = "^(l 5 ) " \pv))- 

The directional coincidence of the p y orbitals with the 
(7 bonds in Fig. 2] identifies them to be of sp type. The 
strength of the a bonds originates from the fact that the 
bands lie 5 to 15 eV below the Fermi energy. 

The physical picture to describe the multicenter bonds 
seems to be much more complicated and it is still under 
investigation (see Appendix IB)) . So far we tried to ana- 
lyze the multicenter bonds using a simple tight binding 



model, which comprises the remaining p x and p z orbitals 
as basis states. But it turned out that this treatment can 
only partially reproduce the conduction bands in Fig. 
probably a larger basis set is needed. 

In Sec. Al we indicated that the puckering has a 
stabilizing effect for the BS. Now we are in a good posi- 
tion to explain this observation: any flattening of the BS 
would cause p x orbitals to interfere with the a bonds and 
eventually destroy them. An analysis of the charge den- 
sity and ELF of a flat BS (see Appendix indeed shows 
that there are no a bonds involved, but only multicenter 
bonds. 

The existence of sp rather than sp 2 hybridization in a 
quasi twodimensional layered structure is somewhat sur- 
prising. Earlier studies of quasiplanar boron clusters^*i& 
still presumed the presence of sp 2 hybridization. We 
think that these assumptions should be reconsidered. 

Finally we want to discuss the two-dimensional Fermi 
surface E F = E BS (k x ,k y ) of the BS in Fig. © It ob- 
viously consists of two contours, which are dispersed 
throughout the Brillouin zone. This clearly shows the 
metallic properties of the BS. 



IV. BORON NANOTUBES 

In Sec. II V Bl we will show that the structure of BNTs is 
strongly related to the structure of the BS, such that the 
latter may be seen as a direct precursor of BNTs. There- 
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TABLE II: Structural data and energies of different isomers of free standing boron nanotubes: (k, I), (n,ra), different 
chiral indices for the same tube type (see Appendix IA lj) : n: number of atoms per unit cell; Isom.: label of isomer; Cy. 
rotational symmetry; ag-_B, ttg^ nal , a grcumferentiai. b oron _b oron bond lengths in axial, diagonal and circumferential direction 
of a nanotube, the superscript a indicates that this bond is a a bond, superscripts a, i and cr, o refer to inner and outer rings, 
respectively; R zb AR: mean radius of a nanotube (Eq. Q) and maximal radial variation (Eq. ©); E 1 ^: cohesive energy of 
a free standing (individual) nanotube (Eq. (Q); -^coiT — ^cSh : this energy is gained when the same nanotube is arranged in a 
bundle (rope). All energies are given in eV/atom and all lengths are given in A. 
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diagonal 
^B — B 
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TTiind 
-^coh 
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(12,0) 
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24 
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0.65 
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+0. 
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(4,4)/(12,0) 


24 
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c 6 




1.69 


1.59", 1.69,1.85 


2.64 
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6.68 
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36 
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1.70,1.74 


1.56", 1.60", 1.71, 1.75 


4.48 
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0.57 


6.74 


+0. 
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1.75 


1.53", 1.76 


4.74 
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0.34 


6.72 






(0,24) 


(8,8)/(24,0) 


48 
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c 6 




1.74,1.75 


1.54"'\1.57"'\1.64"'° 


,1.72,1.74 5.99 
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0.58 


6.81 


+0. 
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fore it will be interesting to try to characterize BNTs 
simply by referring to the BS. The mathematical details 
of such a relation are discussed in Appendix^ and when 
proceeding along these lines, a BNT may be character- 
ized by two numbers (k,l) with fc, I > 0. 

For the ab initio simulation of BNTs we would start 
from a series of initial structures with smooth surfaces, 
which were optimized in a triangular BNT bundle (rope). 
Here the strong tube-tube interactions (see Sec. IIV A|) 
distort the surfaces and naturally induce some puckering. 
The energy of this configuration is ^oh 6 - ^ n OI> der to 
simulate free standing (individual) BNTs we would then 
increase the intertubular distance to 6.4 A, and optimize 
those configurations again while keeping the intertubular 
distances fixed. The energy here is -E^oh- (^coh 6 anc ^ ^coh 
are defined after Eq. (JU)) 48 

All free standing BNTs are shown in Figs. [7[ [51 and 
^2 and the structural data and energies are collected in 
Table ITU Besides their bond lengths and rotational sym- 
metries we also stated the geometrical mean radius of 
each tube R, as well as the maximal radial variation Ai?, 
which were defined as: 



R 
AR 



R n 



R n 



R n 



R — R — R n 



(4) 
(5) 



where R min and R max are the distances of the innermost 
and the outermost atoms from the center of the nan- 
otube, respectively. 



For many (fc, I) BNTs we found more than just one iso- 
mer. Therefore each BNT was also given a Greek index 
which labels different isomers. The latter were ordered 
according to their cohesive energies, i.e., (k,l)a will de- 
note the most stable isomer, (fc, I) (3 would be less stable, 
and so on. 



A. Free standing nanotubes vs nanotube ropes 



In Table [H] the " inter-tubular energy" E r c °^ e 



TTiind 
^coh 



is the energetic difference between a free standing BNT 
and its bundled counterpart. One can see that it varies 
significantly from tube to tube. The intertubular energy 
seems to depend quite strongly on the structure type, 
the relative orientations of adjacent tubes in a rope, and 
the specific type of surface puckering. Furthermore, the 
intertubular distance in different bundles, which was de- 
fined as the minimal separation between two apex atoms 
on adjacent nanotubes, varies between 1.7 and 3.5 A in 
our simulations. 

It is obvious that the tube-tube interaction in BNT 
bundles (ropes) is completely different from what is 
known from carbon nanotubes, where the intertubular 
interaction is of van der Waals type. The latter is cer- 
tainly much weaker, independent of the various structure 
types, and the intertubular distances are always around 
3.4 A. BNTs on the other hand may have covalent inter- 
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FIG. 7: (Color online) The cross sections of different iso- 
mers of a free standing (9,0) zigzag boron nanotube. The 
big spheres stand for the upper atoms and the small ones 
for the lower atoms (with respect to the direction of the tube 
axis) . The a and 7 isomers are the free standing counterparts 
of the (9,0)C and (9,0)B tubes in Ref.E3, respectively. 



tubular bonds ^L2A and this leads to a decent intertubular 
bonding energy that depends quite strongly on structural 
details. 

It is interesting to note that the intertubular energy of 
(0, /) BNTs (armchair types) is significantly higher than 
for (fc, 0) BNTs (zigzag). In Sec. IIVB2l we will try to give 
an explanation for this rather complex bonding scenario. 

At this point, it will be worth noting that the original 
motivation for this paper was a recent study by ourselves, 
where we reported bundled zigzag BNTs that were some- 
what constricted! 7 (we define the concept of constriction 
at the end of Sec. HVB"T]) . We conjectured that this con- 
striction would most likely be caused by the arrangement 
of the tubes in a bundle, where the tube-tube interactions 
will force the tubes to have geometrical shapes different 
from free standing BNTs. Now the free standing counter- 
parts of the constricted (9,0)C and (10,0)C BNTs from 
Ref.|I3are the (9,0)a 49 isomer in Fig. Hand (10,0)/? in 
Fig. [HI To our surprise the constriction would not disap- 
pear after isolating the tube. And even after substan- 
tially deforming the (9,0) a structure by homogeneous 
shrinking, by blowing it up, or by randomly elongating 
atoms out of their equilibrium position with a maximum 
amplitude of 0.3 A, the free standing (9,0)a BNTs would 
always relax to their constricted forms. This finding is 
in clear contrast to our previous hypothesis, and it raises 
the important question where those constrictions finally 
come from. We will try to give an answer to this question 
in Sec. ITvn 
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FIG. 8: (Color online) Cross-sectional view of various isomers 
of free standing (10,0) and (12,0) zigzag boron nanotubes. 
Again the big spheres mark the upper atoms and the small 
ones mark the lower atoms. The (10,0) a and (10,0)/? iso- 
mers are free standing counterparts of the (10,0)B and (10,0)C 
structures in Ref. [ljj respectively. 




FIG. 9: (Color online) Zigzag boron nanotubes and the pres- 
ence of straight a bonds along their axial direction, which 
are indicated by orange (gray) charge density contours at 0.9 
e/A 3 . Due to a lack of stiff a bonds along the circumferential 
direction, this type of nanotube might not be stable. 



B. The structure of free standing boron nanotubes 

1. Zigzag nanotubes 

For zigzag BNTs we found various isomers. Any zigzag 
BNT may be seen as a BS that was rolled up along its 
zigzag direction (see Fig.Elor lT^)) . Thus the linear chains 
of a bonds will lie along its axial direction and they will 
remain straight. These basic bonding properties were 
typical for all zigzag BNT that we studied so far. We 
just show two typical examples in Fig. El Here the bond 
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5 atoms 

FIG. 10: (Color online) Three basic structure elements. The 
cross sections of the most stable zigzag boron nanotubes in 
Figs. □ and |H1 may be composed using these elements only. 



length of the a bonds is quite similar to the bond length 
in the BS, as ag_ B = af£ig = 1.61 A (the only exception 
we found was (9,0)e). 

The tubes (9,0)5, (10,0)7, and (12,0)/? are ideal BNTs, 
which denotes the fact that they were initially con- 
structed by a cut and paste procedure described in Ap- 
pendix an d then reoptimized using ab initio meth- 
ods. Their structure is highly symmetric and we find 
two bond lengths, which are almost identical to the bond 
length in the BS. The puckering height Az = 2AR « 0.8 
is also quite similar to the BS. 

However, an ideal BNT does not seem to be the ground 
state of a real zigzag BNT, and we found less symmetric 
isomers that were higher in cohesive energy. It should 
be noted that zigzag tubes with a smooth surface were 
not considered here because their cohesive energies are 
significantly lower than those of puckered BNTs. As an 
example for the complex shape of zigzag BNTs one may 
study (9,0)e, which is the least stable isomer of all zigzag 
BNTs. (9,0) e has a hexagonal cross section, which prob- 
ably arises from the triangular supercell into which it was 
put. Its sides may be seen as parts of a flat BS, whereas 
the corner pieces are parts of a puckered BS. From Fig.EJ 
we notice that the a bonds along the sides are slightly 
more delocalized than the ones located at the corners. 
This means that any flattening would destabilize the 
sigma bonds, and the whole tube is highly metastable. 
(A similar but squarelike structure was found by Evans 
et a/., 18 which they labeled = (6,6), but we think 
that this structure is highly metastable as well.) Thus the 
question will no longer be if zigzag BNTs are puckered, 
but how they are puckered. 

The cross sections of the isomers with a high cohe- 
sive energy may be built from three basic structure ele- 
ments that are shown in Fig. 1101 The three-atomic struc- 
ture element is directly related to the puckering of the 
BS (compare Fig. |3J) whereas the four- and five-atomic 
elements are just special combinations of three- atomic 
structure elements. We see that the structure of zigzag 
BNTs is strongly related to the local structure of a simply 
puckered BS, but their general cross-sectional geometries 
seem to be more complicated and less symmetric than in 
the simple case of an ideal BNT. This loss of symmetry 
can also be extracted from the spectrum of diagonal bond 
lengths, which are associated with multicenter bonds. 



Those bond lengths are not equal to a^^ nal of the BS 
(see Table H|}, but span a whole range a^^ nal « 1.7 — 1.9 

A. 

Some of the most interesting structures are (9,0)a, 
(9,0)/?, and (10,0)a, which have cross sections that are 
far from being circular. Nonetheless they exhibit high co- 
hesive energies. Because of the observed unusual shapes 
of zigzag BNTs we assume that the multicenter bonds 
obviously possess a high directional flexibility, but at the 
same time they are also very stiff (C x = 0.42 TPa in Ta- 
ble [I]) . Therefore it seems as if these bonds have some 
jointlike properties, i.e., they are easy to turn, but hard 
to tear. 

In the following we will call a zigzag BNT constricted, 
if it is composed of several five- atomic structure elements. 
In our work the (9,0)a and the (10,0)/? isomers are con- 
stricted. A constricted zigzag BNT was also found by 
Evans et ali& There it is labeled as a (i, j) = (8, 8) nan- 
otube, and it corresponds to our (10,0)/? structure with- 
out the two horizontal three-atomic elements. 



2. Armchair nanotubes 

When rolling up a BS along its armchair direction, the 
puckered sheet (see Fig. 0J will be transformed into a 
tube that has inner and outer rings, and the a bonds 
will lie along its circumferential direction. On the outer 
rings the length of the a bonds will be increased and on 
the inner rings their length will be reduced. In Fig. 
we see that for three systems discussed in this study, the 
o~ bonds do really lie along the circumferential direction, 
and for the (0,18) and (0,24) systems an inner and an 
outer ring can clearly be identified. 

In contrast to zigzag BNTs, for the armchair types 
we did not find several isomers, and we just discuss one 
ideal BNT, which is the (0,18)/? isomer. In analogy to 
zigzag BNTs, we found that this ideal BNT corresponds 
to a local energy minimum, and the (0,1 8) a isomer of 
lower symmetry is 0.02 eV/atom more stable. 50 The lat- 
ter has o~ bonds solely along the inner rings, where the 
bond lengths are 1.56 and 1.60 A. Along the outer ring, 
where the B-B distances (1.71 and 1.75 A) are signifi- 
cantly longer, the curvature effect has destroyed the a 
bonds™ 

The (0,24) system has similar properties, but here the 
curvature is smaller, and there are six additional weak a 
bonds along the outer rings with a bond length of 1.64 A. 
For even larger radii we expect the outer rings of armchair 
BNTs to develop a bonds between every single atom. 

The radius of the (0,12) BNT is quite small, which 
makes it extremely difficult for the structure to align its 
o~ bonds. We see that this tube possesses a different ge- 
ometry, and even along the stiffer rings there are six in- 
stead of 12 a bonds. It is obvious that for armchair BNTs 
with smaller and smaller radii, the curvature effect will 
successively destroy the circumferential a bonds. For the 
smallest possible BNTs there will probably be no a bonds 
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(0,12) a (0,1 8) a (0,1 8) p (0,24) a 



FIG. 11: (Color online) Top and side view of various free standing armchair boron nanotubes and the presence of a bonds, 
which are indicated by orange (gray) charge density contours at 0.95 e/A 3 . All armchair nanotubes have bent a bonds along 
the circumferential direction, which basically generate the strain energies of the tubes. The black bar on the right indicates 
the height of a supercell in axial direction that was used for our simulations; for aesthetical reasons we actually displayed three 
identical units cells. 



at all, and the surface of the tube will become smooth. 
This agrees with earlier studies by ourselves*^ and with 
the work of Evans et al., 18 where some armchair BNTs 
of small radii were studied and found to be smooth. 

Any destruction of circumferential a bonds within 
armchair BNTs of small radii will release electrons that 
can alter their chemical properties. In Sec. II V Al we ob- 
served that the intertubular energy for armchair BNT 
ropes is much higher than for zigzag BNT ropes. Now a 
possible explanation would be that the released electrons 
in armchair BNTs induce an enhanced reactivity. In a 
rope of BNTs, this enhanced reactivity will lead to strong 
intertubular bonding for armchair BNTs of small radii. 
In zigzag BNTs the reactivity is lower, as a maximum 
number of a bonds can always be achieved, due to the 
fact that curvature effects will not be able to weaken the 
axial (7 bonds. Therefore we hypothesize that small sized 
armchair BNTs will have a higher reactivity than zigzag 
BNTs, and that this reactivity will further decrease with 
increasing radii. 

This reactivity, which leads to the formation of inter- 
tubular bonds in BNT ropes, could be very useful when 
trying to embed BNTs into polymers, 18 where strong 
chemical bonds between the nanotubes and the polymer 
matrix are needed in order to improve the mechanical 
properties of the composite. 



C. Strain energy 

Let us now compare the cohesive energy of every BNT 
(£coh f rom Table|n]) with the cohesive energy of the puck- 
ered BS (i^h from Table [I]). This energy difference will 
be called strain energy: 

E strain (k,l) = Efi-EZi(k,l). (6) 

It is the amount of energy that is needed to roll up a BS 
into a BNT. The microscopic origin of the strain energy 
in nanotubes are bent a bonds along the circumferen- 
tial direction of the tubes. These bonds have a strong 
tendency to jump back into a straight orientation, which 
generates a tension that may be quantified by the strain 
energy of the systems. Such a tension will stabilize the 
tubular shape, or to put it more clearly: it will make the 
nanotube round. 

The strain energies of different (fc, I) BNTs as a func- 
tion of their mean radii (Eq. Q) is plotted in Fig.^J For 
the sake of comparison we also show the universal strain 
energy curve for carbon nanotubes. We call it universal 
because the strain energy only depends on the radius, 
but not on the chiral angle (chirality) of the nanotubes: 

^strain = ^strain (^)* 

As the BNTs are all puckered, there is some variability 
in the proper choice of a mean tubular radius. Since the 
strain energy is related to the position of the a bonds, it 
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FIG. 12: (Color online) The strain energy of a isomers as a 
function of the mean radius R (Eq. ©); for armchair boron 
nanotubes we used Ra (Eq. Q). In orange (gray) we show 
the universal strain energy curve for carbon nanotubes (□); 
the energy obviously depends on their radii, but not on their 
chiral angles. For armchair boron nanotubes (O) we find a 
similar curve, but those boron tubes have more strain energy. 
For zigzag boron nanotubes (A) we cannot really plot a strain 
energy curve, as different nanotubes of different radii are al- 
most isoenergetic. Ideal zigzag boron nanotubes (O) have 
less strain energy than their armchair counterparts, but they 
are met ast able. 



makes sense to define the mean radius of armchair BNTs 

as 



R 



(7) 



Here R™ m and i?™ ax are the distances of the innermost 
and the outermost atoms sharing a bonds, which is mea- 
sured from the center of the nanotube. 

In earlier works we studied the elasticity of armchair 
BNTs with a tight-binding method 13 and reported a typ- 
ical strain energy curve lying below the one of carbon 
nanotubes H Now, using an ah initio method, we also 
found that armchair BNTs have strain energy, but it is 
higher than for carbon nanotubes. 

Different ideal zigzag BNTs in Fig. [El have rather low 
strain energies. Here none of the a bonds has to be bent, 
and the strain energy should only come from the multi- 
center bonds. But those ideal BNTs are metastable, and 
isomers of lower symmetry have higher cohesive energies. 
Thus for the zigzag a isomers no strain energy curve may 
be plotted as they are more or less isoenergetic. It seems 
that zigzag BNTs can release some or all of their strain 
energy by lowering their symmetry and undergo internal 
deformations (see also Ref. 18), possibly mediated by the 
jointlike properties of the multicenter bonds. 

In summary we see that the strain energy in BNTs is 
mainly caused by bent a bonds lying entirely (armchair) 
or only partially (chiral BNTs) along the circumferential 
direction. The multicenter bonds are always present, but 
they seem to have no serious effect on this. The apparent 



absence of strain energy in zigzag BNTs is caused by the 
fact that the linear a bonds lie along the axial direction, 
only. But without smoothing bonding strains, the zigzag 
tubes are free to take a multitude of cross-sectional mor- 
phologies. This explains the number of different isomers 
that we found for (9,0), (10,0), and (12,0) zigzag BNTs 
and their bizarre shapes. The constriction of zigzag BNT, 
first reported in Ref. [I?! is a clear consequence of the ab- 
sence of strained bonds within zigzag BNTs. Armchair 
BNTs in turn, which are geometrically stabilized by their 
strain energy, do not seem to have this kind of isomerism. 

Chiral BNTs may be thought of as a certain combi- 
nation of structural elements from armchair and zigzag 
tubes defined by a certain chiral angle. Therefore we 
suppose that there will be a separate strain energy curve 
for every chiral angle lying in-between the armchair and 
the zigzag curves. The strain energies themselves will 
depend on the radii and on the chiral angle of a BNT: 
^strain = ^strain C^j^)- This seems to be a unique prop- 
erty among all nanotubular materials reported so far. 

But it remains open whether the strain energy of 
zigzag BNT will be completely absent, or just signifi- 
cantly smaller than for armchair BNTs. The present re- 
sults are in favor of the former hypothesis. As carbon 
nanotubes with large diameters (and very small strain 
energies) are susceptible to a structural collapse j^lLSi it 
is possible that without a significant amount of strain 
energy the zigzag nanotubes could be geometrically un- 
stable. Given some thermal fluctuations or strain they 
might collapse just like big diameter carbon nanotubes. 
However, such a collapse might also be prevented by a 
possible energy barrier, which should be proven to be 
absent in order to support this collapse hypothesis. 

Finally we want to point out that the constriction of 
zigzag BNTs could be an important intermediate mecha- 
nism during the possible collapse of a zigzag BNT, which 
might allow for the formation of B12 icosahedra, that 
are the basic building blocks of all previously known bulk 
boron structures. The five-atomic element (see Fig. ITU)) 
forms part of an imaginary zigzag 6-ring, similar to the 
six apex atoms of a B12 icosahedron, as seen along each 
of its threefold axestf 



V. SUMMARY AND CONCLUSIONS 

In this paper we studied a number of different struc- 
ture models for broad boron sheets (BSs). All of them 
are metallic, and we found that for a 16 atom supercell, 
the model with a simple up-and-down puckering will be 
the most stable one. Large quasiplanar boron clusters 
with a similar structure (B22^ B48 and Bgg 16 ) were 
already reported before. They may now be understood 
as a first indication for the onset of periodicity in finite 
layered boron systems, and thus they are an independent 
confirmation of the current structure model. 

A fiat BS has a rather high stiffness, and it seems to 
be held together primarily by multicenter bonds (see Ap- 
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pendix lB|) . Although the sheet is less stable than previ- 
ously known bulk phases of boron, as shown here and 
elsewhere ^2ii& the model sheet could be the ideal theo- 
retical tool for studying complex multicenter bonds. 

After describing the lattice structure of the stable BS, 
we have analyzed its band structure, the corresponding 
charge densities, and the electron localization function. 
This would lead to the following preliminary picture of 
the chemical bonding: on the one hand the sheet is held 
together by homogeneous multicenter bonds, on the other 
hand there are linear sp hybridized a bonds exclusively 
lying along the armchair direction of the sheet. The exis- 
tence of sp hybridization in quasiplanar BS is somewhat 
surprising given the fact that earlier studies would always 
claim sp 2 hybridization. The rather anisotropic bond 
properties of the sheets lead to different elastic moduli 
C x and C y for stretching the BS in the x and in the y di- 
rection. Furthermore puckering of the BS, which breaks 
the triangular symmetry, may be understood as a key 
mechanism to stabilize the sp a bonds. Our results in- 
dicate that the sheet analyzed in this study is the boron 
analog of a single graphene sheet, a possible precursor of 
boron nanotubes (BNTs), and we wonder whether broad 
BSs might exist in nature. 

Constructing BNTs from the BSs by a cut and paste 
procedure will generate ideal BNTs (see Appendix EJ . 
Because the underlying two-dimensional lattice structure 
is rectangular rather than triangular or hexagonal, it fol- 
lows that the chiral angle ranges from 0° to 90° (0 = 0°: 
zigzag, 6 = 90°: armchair), and that chiral BNTs do not 
have an axial translational symmetry. We therefore pre- 
dict the existence of helical currents in ideal chiral BNTs 
( Appendix IA 2|) . Furthermore we presented a band the- 
ory for ideal BNTs, employing their helical symmetry, 
and showed that all ideal BNTs are metallic, irrespective 
of their radius and chiral angle (Appendix IA 3|) . BNTs 
could therefore be perfect nanowires, superior to carbon 
nanotubes. 

In an independent study of armchair and zigzag BNTs 
we found that ideal BNTs do not form the ground state 
of BNTs, and we identified structures of lower symmetry, 
which are higher in cohesive energy. The symmetries of 
real BNTs still remain to be determined, and the ideal 
BNTs may be seen as rather close approximants to real 
BNTs. 

We also found that all BNTs, except small radius arm- 
chair types, had puckered surfaces, and a bonds along the 
armchair direction of the primitive lattice. The existence 
and mutual orientation of these a bonds is crucial to un- 
derstand the basic mechanical and energetic properties of 
BNTs because the strain energy of the tube is mainly gen- 
erated by bending those a bonds. The multicenter bonds 
seem to have no real effect on the strain energy. They are 
likely to have joint like properties (they are easy to turn 
but hard to tear), which allows for a certain flexibility of 
these bonds, and any bonding strain could immediately 
be released through internal relaxations™ 

We showed that armchair BNTs, where the a bonds lie 



along the circumferential direction, will have rather high 
strain energies, whereas zigzag BNTs, where the a bonds 
will lie along their axial directions, will have nearly van- 
ishing strain energies. Thus BNTs have a strain energy 
that depends on the nanotube's radius R as well as on the 
chiral angle 0: Ef train = Ef tTahL (R, 0) . We suppose that 
there will be an individual strain energy curve for every 
chiral angle lying between the armchair and the zigzag 
curves. This is a unique property among all nanotubular 
materials reported so far. 

This intriguing feature could even allow for some struc- 
ture control in nanotechnology. For carbon nanotubes, 
the strain energies do not depend on their chiralities 
(E^ train = E^ ra[n (R)), and thus the experimentalists may 
control the radius, but not the chirality of carbon nan- 
otubes, although the latter will determine the electronic 
properties of such materials. With the experimental tech- 
niques at hand today one might be able to walk along 
the energy axis by tuning the reaction conditions, and 
along the radius axis by synthesizing nanotubes within 
porous templates with well defined pores sizes. This way 
it could be possible to synthesize BNTs of a specific type 
only. The connection to carbon nanotubes may occur 
via intramolecular junctions, 37 allowing for a controlled 
layout of carbon nanotubes as well. 

The rather low strain energies in zigzag BNTs lead to a 
whole bunch of possible structural isomers, as a nanotube 
without any significant amount of strain energy will not 
be able to maintain a circular cross section. This can lead 
to a certain constriction of zigzag BNTs^i and we even 
hypothesize that zigzag BNTs could be too unstable to 
really exist out in nature, provided that there will be no 
significant energy barrier left to prevent a collapse. 

Armchair BNTs on the other hand are geometrically 
stabilized by their strain energies, but for armchair BNTs 
of rather small radii, the BNTs are unable to maintain 
a puckered structure necessary to align the circumfer- 
ential a bonds. In agreement with earlier studies 5 ^^ 
we expect them to flaten out and build up a smooth 
surface. Furthermore, we hypothesize an enhanced re- 
activity of small radius armchair BNTs in comparison to 
zigzag BNTs, which could be useful for embedding BNTs 
into polymers. 18 
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FIG. 13: (Color online) The triangular (t), the rectangular 
(r), and the honeycomb- derived (h) primitive cells that are 
used to characterize boron nanotubes. They contain one, two, 
and three atoms, respectively. Only the rectangular cell may 
properly describe the puckering of the boron sheet (indicated 
by black and gray atoms in the background). 



APPENDIX A: THE MATHEMATICAL 
DESCRIPTION OF IDEAL BORON NANOTUBES 

1. Wrapping vector 

The geometrical construction of BNTs from BSs is 
similar to the construction of carbon nanotubes from a 
graphene sheets the basic tubular structure is charac- 
terized by a wrapping vector W that defines a rectangu- 
lar area on the BS, which is rolled up to a cylinder such 
that W becomes the circumference of the nanotube and 
its radius will be R = \ W\/2tt (see Fig. lHJ) . We will call 
any BNT, whose structure may be described by such a 
construction, an ideal boron nanotube. 

Due to the fact that a proper structure model for BS 
was missing for a long time, there remains some confusion 
in the literature about a proper reference lattice struc- 
ture. In the work of Cabria et alm& and in earlier works 
by usiiii (and in full analogy to the construction of car- 
bon nanotubes) the BNTs are related to a honeycomb 
lattice and denned the wrapping vector W h as 

W h = (n, m) = na\ + ma\, (Al) 

a\ 2 are the primitive vectors of a honeycomb lattice and 
n, m are integers. Here each unit cell has one additional 
atom at the center of the honeycombs, thus consisting of 
three rather than two atoms (see Fig. IT3|) . Gindulyte et 
al.fi Evans et al.ML and some earlier work of ours 13 relate 
their BNTs to the simple triangular lattice, having only 
one atom per unit cell: 

W t = (i,j)=ia\+ja t 2 , (A2) 

a\ 2 are the primitive vectors of a triangular lattice, and 
z, j are integers. W h and W l can be transformed into 
each other by using 53 

(n, m) i — ^ (i, j) = (n + 2m, n — m), (A3) 
h-> (n,ra) = + 2j, i-j). (A4) 

From Fig. EH we see that both definitions are based 
on primitive vectors, which have different orientations. 54 




zigzag 



FIG. 14: (Color online) The geometrical construction of an 
ideal boron nanotube from a boron sheet: the red (gray) area 
is cut and rolled up such that W r will become the circumfer- 
ence of the nanotube. O is the origin, W r is the wrapping 
vector, T is the translational vector, 6 is the chiral angle 
measured with respect to the zigzag direction, ai,2 are the 
primitive vectors of the underlying rectangular lattice, and A 
and B are the lattice constants (see text). The puckering of 
the boron sheet is indicated by black and gray atoms in the 
background. The zigzag and the armchair directions are per- 
pendicular to each other. This figure corresponds to W r = 
(5,3) and A/B = >/3, which implies T = (-1, 5). 



This leads to the rather unsatisfactory situation that 
armchair and zigzag directions are swapped in both de- 
scriptions (see Table IIII for example). Cabria et al. found 
that all (n, 0) zigzag and all (2n, 2n) armchair BNTs have 
puckered surfaces, while the (2n + l,2n + 1) armchair 
tubes shall be smooth due to the fact that an odd num- 
ber of boron rows along the tube surfaces does not allow 
for the formation of the simple up and down puckering. 19 
We think that these results are not an intrinsic property 
of BNTs, but rather a consequence of an unsuitable ref- 
erence lattice system that is unable to properly describe 
the puckering of the boron sheet, see Fig. [E3 Further- 
more, the puckering will break the hexagonal symmetry 
underlying the honeycomb and the triangular lattice. 

Therefore we convinced ourselves that these descrip- 
tions are not really appropriate to classify BNTs. On 
the basis of the current BS model we would like to put 
forward a different way of describing BNTs, based on a 
rectangular lattice underlying the two-dimensional struc- 
ture of the BS. 

We define the wrapping vector W r as 

W T = (k,l) = ka\+la\, (A5) 

fc,Z are integers, and a\ = A(1,0) and a\ = 5(0,1) 
are the primitive vectors of the rectangular lattice (see 
Figs. ESI and EI); A and B are the lattice constants from 
Table [I] In analogy to the Dresselhaus construction for 
carbon nanotubes we define the chiral angle as the 
angle between the vectors W r and a\, i.e., is measured 
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with respect to the zigzag direction coinciding with a\ 
(see Fig. HI)) . 

The categorization of BNTs will be different from other 
classification schemes because the reduced symmetry of a 
BS increases the number of possible types of nanotubes, 
as the range for the chiral angle will be 0° < 6 < 90°, 
and for the chiral indices (k,l) we find that k,l > 0. 
Zigzag BNTs will now correspond to = 0° and (/c, /) = 
(&,0), and armchair BNTs will correspond to = 90° 
and (fc,Z) = (0,/). 

W h and W t cannot directly be converted to W r , as 
they are defined for lattices with different symmetries. 
For the achiral types, on may use the following analogy 
(for examples see Table ITT)) 

zigzag: (/c,0) r <-> (k,0) h (A6) 
armchair: (0, l) r (l/3,l/3) h 

- (hoy. 

2. Translational vector 

The tubular unit cell of an ideal BNT, being the red 
(gray) area in Fig. El m &y be defined properly by a wrap- 
ping vector W T (Eq. J23J)) and the so-called translational 
vector T, which is perpendicular to W r : 



T = (t u t 2 ) 
t 2 



t 1 a r 1 +t 2 a r 2 , (A7) 

J — numerator (/) : k ^ 

\ 1 : k=0 

J denominator(/) : k ^ 

\ : k=0 

^duce (^) . 



ti, t 2 are integers and reduce (r) should indicate that the 
fraction r must be reduced before determining its numer- 
ator and denominator. 

Let us consider the length of the translational vector 
T. For the achiral BNTs |T| is particularly small: for 
all (fc,0) zigzag types we have T = (0,1), and for (0,/) 
armchair BNTs T = (1,0). For the chiral types T de- 
pends on the ratio B 2 /A 2 (see the last line of Eq. (S3)- 
Using A = 2.819 and B = 1.602 we obtain reduce 
(B 2 /A 2 ) = 2566404/7946761. Therefore the coefficients 
t\ and t 2 are really huge numbers, which means that |T| 
becomes macroscopically large. For A and B chosen as 
above, the estimated length of T for all chiral BNTs will 
be in the mm range. Imposing some additional symmetry 
constraints by relating the lattice constants will imme- 
diately remedy this problem. For example by choosing 
A = V3B, fraction(£ 2 /A 2 ) = 1/3, i.e., |T| will be re- 
duced to just a few lattice constants (this case was used 
to generate Figs. HH and 113)) . So for the chiral BNTs the 
specific ratio B 2 /A 2 determines the length of the trans- 
lational vector. 




FIG. 15: (Color online) Two different ways of "building up" a 
nanotube: the tubular unit cell in light gray (see also Fig.[2J 
is repeated along the nanotube's axis, which lies parallel to 
T. The helical unit cell in red (dark gray) is translated along 
spirals (represented by the dotted lines) on the surface of the 
nanotube; it is defined by the helical vector H and vector 
K . It holds W v = K + L. Here H (0,1), K (2,0) 
and L = (0,9), and therefore W r = (2,9). The length of 
T — (—3, 2) was artificially reduced by choosing A/B — y/3. 



Boron compounds usually have a whole set of differ- 
ent B-B bond lengths, which means that boron does not 
necessarily favor highly symmetric arrangements. The 
bond lengths are more flexible than for typical covalent 
elements like carbon, and the lattice constants A and 
B of the BS cannot really be seen as fixed parameters; 
they will have slightly different values in BNTs. Further- 
more, the broken planar triangular symmetry of the BS 
is rather typical for boron, and we should expect that for 
ideal chiral BNTs, even with different values of A and 5, 
the translational vector might still be large. 

To summarize: any departure from triangular sym- 
metry in the BS will create chiral BNTs, which have 
macroscopically large translational vectors, and achi- 
ral types, where |T| is of the order of the lattice con- 
stants. Thus achiral BNTs (armchair and zigzag) have a 
one-dimensional translational symmetry along the tube's 
axis, which is not present in chiral BNTs. For the latter 
it might be better to think in terms of helical (chiral) 
symmetries only. Therefore we predict the existence of 
helical currents in ideal chiral BNTs. Such currents could 
lead to very interesting physical effects such as strong 
magnetic fields 4 - and self-inductance effects leading to 
an inductive reactance 41 of chiral BNTs. 



3. Band structure 

Within the limit of large nanotube radii, where curva- 
ture effects are small, one may derive the one-dimensional 
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band structure of an ideal BNT E^(k') by a zone- 
folding technique^ starting from the two-dimensional 
band structure of a BS E BS (k x , k y ). Given the absence 
of translational symmetry in ideal chiral BNTs, we have 
to base our zone-folding theory on the helical symmetry 
of BNTsiii4 

Figure fTTTI illustrates that besides constructing a BNT 
by repeating a tubular unit cell one can also build a nano- 
tube by repeating a helical unit cell along a spiral wind- 
ing around the surface of the tube. The direction of this 
spiral is given by the helical vector H 42A3 (in Ref. 39 it 
is called the symmetry vector which, when uncoiled 
into a plane, defines the direction of a translational sym- 
metry (see Eq. (|A8j) and thereafter). The helical unit cell 
is specified by H and the vector K JL H. Furthermore, 
we define a vector L \\ H, such that W v = (fc, I) = K + L 
(see Fig. IT5J) . 

The helical wave functions are restricted by the follow- 
ing criteria: 

¥ Mfc /(r + Jr) = ^^(r)exp(ik'\H\), (A8) 
V^'ir + W 1 ) = *^/(r). (A9) 

Equation (|A8|) defines a one-dimensional Bloch state with 
—tt/\H\ < k' < tt/\H\ and imposes the condition that k' 
has to be parallel to the reciprocal vector of H. Equa- 
tion (|A9|) is the tubular boundary condition. In order 
to construct the helical wave functions we have to 
use the wave functions of the BS \I/^ s (r) which have the 
Bloch property: 

* BS (r + R) = exp(*fc • iZ)^ s (r), (AlO) 

where R is a vector of the Bravais lattice formed by a\ 
and a\. Since the vectors H and W r are also elements 
of such a Bravais lattice, Eq. (|A8|) will automatically be 
satisfied, and Eq. (|A9|) together with Eq. (|A1Q|) will yield 

l = exp[i(fe. W v )\. (All) 

In order to proceed, we now have to define the direction 
of H, which may be any Bravais lattice vector. 55 By 
choosing H = T we recover the case of a tubular unit cell, 
as described above and in Ref. But in order to make 
the calculation as simple as possible we assign H = a\ = 
(0,1). Then it follows that K = (fc,0) and L = (0, 1) 
(see Fig. lTSJ) . As H || y we have to choose k' = k y . After 
inserting Eq. (|A11|) into E BS (k x , k y ) we finally obtain the 
zone-folded band structure of ideal (fc, I) (k ^ 0) BNTs 
as 

E i kJ) w = EBS (u»-u k '> k ')> (A12) 

~B <k <B> 
fi = 0, • • • , k — 1 

Equation (|A12|) will break down for (0,1) armchair 
BNTs, due to a chiral index k = 0. But as mentioned 
before, we are free to choose the direction of H, and in 




FIG. 16: (Color online) Properties of a flat boron sheet: (a) 
The two-dimensional band structure, (b) Black lines indicate 
the triangular unit cells, black spheres are boron atoms, and 
the orange (gray) contours show the electron localization func- 
tion (ELF) at contours of 0.7. We observe a simple network 
of two- and three-center bonds. 

such a case we use H = a\ = (1,0) and have k' = k x . 
We thus obtain 

E^\k')=E* s (k',^), (A13) 

7 < k < -J' 

A A 

11 = 0, 1 

Unfortunately we do not have an analytical band struc- 
ture of the BS E BS (k x , k y ), yet. But to decide whether a 
certain ideal BNT is metallic or not we can simply zone- 
fold the BS's Fermi surface given in Fig. El We did so 
and found that all ideal BNTs are indeed metallic, irre- 
spective of their radius and chiral angle. The only ideal 
BNTs that are not metallic are the (0,1) and the (0,2) 
types. But these structures are highly unrealistic and we 
can safely rule them out, as they are not even covered by 
the Aufbau principle. 3 

APPENDIX B: A FLAT BORON SHEET 

The lattice structure, the cohesive energy, and the elas- 
tic moduli of the flat boron sheet - model (a) - can be 
found in Table |I] The elastic modulus of Cn ~ 750 GPa 
is comparable to graphite. The electronic charge density 
is nearly uniform in the interstitial region, and the band 
structure (see Fig. lTCT a)) is similar to the band structure 
of a free electron gas. These results seem to indicate some 
metallic bonding, as pointed out by Evans et al.^ but 
such a picture cannot really account for the planarity and 
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FIG. 17: (Color online) A kinked boron sheet based on a 
supercell (thin lines) that contains 16 atoms (see text). Apart 
from the kink its surface is slightly puckered. 

the high elastic modulus of the flat BS. A different qual- 
itative picture of the chemical bonding is obtained after 
looking at the electron localization function^ (ELF) in 
Fig. rTrjT b). Here we observe a simple network of two- and 
three-center bonds being less localized (ELF w 0.7) than 
typical sigma bonds (ELF « 0.9), which are absent here. 
Thus the flat BS seems to be held together predominantly 
by multicenter bonds similar to the ones found in pure 
boron compounds. The chemical understanding of these 



multicenter bonds is still very limited. We think that, de- 
spite of its apparent met ast ability, model (a) could be an 
ideal theoretical tool to extend our present understand- 
ing of the nature of multicenter bonding in boron. 



APPENDIX C: A KINKED BORON SHEET 

In Sec. IIII Al we described the optimization of randomly 
puckered BSs. It is surprising that despite the high com- 
plexity of the boron energy landscape for small boron 
clusters, which are known to have many local minima, 
these runs seem to have only two possible " attractors" . 
One is model (b) - a simply puckered BS - the other 
is the kinked BS displayed in Fig. [T7| The kinked BS 
has a metallic density of states and a cohesive energy of 
6.86 eV/atom, which is somewhat intermediate between 
model (a) and model (b). We think that this structure 
is likely to be an artifact of the finite size of the super- 
cell that we used for the simulation runs, but being an 
" attractor" of the optimization runs it is still interesting 
enough to be mentioned here. 
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